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Abstract: We revisit BPS solutions to classical N = 2 low energy effective gauge theories. 
It is shown that the BPS equations can be solved in full generality by the introduction of 
a Hesse potential, a symplectic analog of the holomorphic prepotential. We explain how 
for non-spherically symmetric, non-mutually local solutions, the notion of attractor flow 
generalizes to gradient flow with respect to the Hesse potential. Furthermore we show that 
in general there is a non-trivial magnetic complement to this flow equation that is sourced 
by the momentum current in the solution. 
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1. Introduction and summary 

Theories with N = 2 supersymmetry in 4 dimensions seem to strike a unique balance 
between analytical control and non-trivial physics. This is most famously illustrated by 
the work of Seiberg and Witten [1, 2] in gauge theory and that of Maldacena, Strominger 
and Witten [3] in the presence of gravity. What makes N = 2 theories so special and 
suitable to an exact analysis is the presence of BPS states. Finding the spectrum of stable 
states of a given physical theory is in general a formidable problem, especially at strong 
coupling. The simpler problem of computing the BPS spectrum remains very challenging, 
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but recently a number of new insights have led to tremendous progress that suggests a 
complete and exact solution to this and other problems in N = 2 theories might be within 
reach. See e.g. [4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14] for a small selection out of the everyday 
growing literature. 

BPS states have been studied through a wide variety of methods. Among others one 
can think of them as non-Abelian dyons [15], string- webs [16], WKB curves [7], attrac- 
tors [17, 18] or ground states in (quiver) quantum mechanics [19, 13, 14]. Each of these 
representations has its merits and limitations. 

A particularly interesting regime is that of low energy. Studying BPS states in this 
limit has many advantages. Through electromagnetic duality we have an exact description 
at each value of the coupling [1]. Furthermore all gauge interactions are Abelian and 
this simplifies the theory a lot. Finally this description is directly formulated in physical 
4d space-time, making it very intuitive. The formalism has mainly been developed in 
supergravity [18, 20, 4], since there the UV (string) theory is only partially understood. 
This study of BPS solutions to M = 2 supergravity led to an interesting conjecture. Denef 
proposed that BPS states and their stability are encoded in so called split attractor flows [18, 
4]. To prove this conjecture we would need a much better and more detailed understanding 
of Type II string theory compactified on a Calabi-Yau manifold. 

In the case of N = 2 gauge theories the UV theory is under much better control and 
in certain cases the BPS spectrum is fully understood, in one or more of the descriptions 
mentioned above. This suggests that an analog of the split attractor flow conjecture in 
gauge theory could be subjected to precise tests. This could furthermore lead to a precise 
formulation and possibly a proof, or if not, at least an understanding of its limitations. 
This paper can be seen as a first step towards formulating split attractor flows in N = 2 
gauge theory. 

BPS solutions to low energy M = 2 gauge theory have been considered in the literature 
before [21, 22, 23, 24, 25]. However, none of these have made a completely general analysis. 
Either one concentrated on some specific (simple) theories or the focus was on spherically 
symmetric solutions to an arbitrary N = 2 theory. This last assumption is especially strin- 
gent, as it is well known that much of the spectrum manifests itself as solutions that consist 
of multiple centers carrying mutually non-local charges and breaking spherical symmetry 
[26, 18]. 

In this paper we present the general solution to the BPS equations of any low energy 
J\f = 2 gauge theory without matter. In section 2 we will review some essential properties 
of J\f = 2 gauge theories and their low energy description in terms of r massless vector 
multiplets and the prepotential T . We then write down the BPS equations in section 3, 
and show how a general solution is determined in terms of 2r real harmonic functions 1 . 
The relation to the original r vector multiplets is shown to involve a Legendre transform 
of the imaginary part of the prepotential that is known in the supergravity literature as 
the Hesse potential S [27, 28, 29]. 

1 The precise description of the sources for these harmonic function and the interactions among them 
is an interesting problem that we shortly discuss in section 4.3 but of which we leave a general detailed 
treatment to future work. 
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This Hesse potential first appeared in discussions of the special geometry of the theory 
[30, 31, 32, 33, 34, 35]. It takes its name from the fact that in Darboux coordinates the 
metric on moduli space is nothing but the Hessian matrix of S. We review these facts 
in the first part of section 4 and appendix A. Since both the Hesse potential and the 
Darboux coordinates Y a give a duality covariant description of the special geometry we 
can use them to give a duality covariant description of the complete theory, including the 
U(l) gaugefields. We formulate the duality covariant Lagrangian (4.10) and its variation 
principle in the second part of section 4. In appendix B we show in detail the equivalence 
between this formulation and the standard one. In this section 4 we also rederive the BPS 
equations in a manifestly duality invariant form and at the end comment on the description 
of sources and the constraints on their positions. 

In section 5 we show how the attractor equations for the central charge, well estab- 
lished for spherically symmetric solutions, generalize to gradient flow equations determined 
by a central charge vector field. It is here that the Hesse potential finds its physical in- 
terpretation: evaluated on the solution it is the potential for this gradient flow, that we 
therefore refer to as Hesse flow (5.8). Finally we show how the BPS equations imply a 
rather remarkable set of equations for the central charge vector field. If we interpret its 
real part as an 'electric' field and its imaginary part as a 'magnetic' field then it satisfies 
Maxwell's equations (5.14, 5.15). What is particularly intriguing is that the equations are 
sourced by the energy- momentum 4-current. In particular these equations imply that the 
Hesse potential evaluated on the solution equals the gravitational potential of that solution 
(i.e. it solves the Poisson equation sourced by the energy density). Furthermore in cases 
without spherical symmetry and mutually non-local charges there is an additional magnetic 
equation sourced by the momentum current present in the solution. 

2. M = 2 gauge theory at low energy 

In this section we shortly review some notions in M = 2 gauge theory and in particular 
their low energy description. We will only consider results that are relevant to the rest of 
the paper and use this section mainly to fix our notation. 

In this paper we will restrict ourselves to the vector multiplet sector of N = 2 rigid 
supersymmetry in 4 dimensions. We expect our results to remain essentially unchanged, 
at least qualitatively, when hypermultiplets and/or the gravity multiplet are included and 
some crucial notions are appropriately adjusted to those cases. 

2.1 UV origin 

The N = 2 vector multiplet consists of a gauge field A^, a complex scalar (f> and a Dirac 
spinor tp, all valued in the adjoint representation of some gauge group G. The extended 
supersymmetry forbids the presence of a superpotential and so the scalar potential arises 
purely from integrating out auxiliary superfields and is given by 
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In a generic vacuum the gauge group G will be spontaneously broken to its maximal torus 
U(l) r , with r = rank G. Furthermore, since the flat directions of the potential (2.1) are 
in one to one correspondence with generators of the Cartan subalgebra of G, there are 
also precisely r massless scalar fields in such a vacuum. This implies that at low enough 
energies the theory is effectively described by the dynamics of r massless U(l) vector 
multiplets (A A , <f> A , i(j a ). 

J\f = 2 supersymmetry is restrictive enough to fully determine the two-derivative low 
energy effective action up to a single holomorphic function, the prepotential T. In M = 2 
superspace the action takes the simple form 



S = ImJ d A xd A 6F(W). (2.2) 



The precise form of F depends on the details of the UV theory and encodes the effect 
of integrating out the various massive fields. Using arguments based on electromagnetic 
duality, Seiberg and Witten [1] provided a geometric way to efficiently and exactly compute 
F in the case G=SU(2). This result has since been generalized to include a wide variety 
of gauge groups. Some of the state of the art can be found in [36], together with a nice 
overview of different techniques and further references to the literature. 

In this paper we will make no further assumptions about F and start directly from the 
low energy effective action, leaving the rank r and prepotential F as arbitrary input param- 
eters. Readers interested in more concrete applications can directly apply our formalism 
to any specific r and F obtained from a given UV theory. 

2.2 Low energy bosonic Lagrangian 

Let us now take a closer look at these low energy N = 2 effective actions. The bosonic 
action following from (2.2) is 



S = - f d 4 xlm TAB (±F*(F B ^ + i*F B »») + ^<f> A d^ B ) 



(2.3) 



As explained above, the fields in this Lagrangian are r complex scalars 4> A and r U(l) 
gaugefields A A , labeled by an index A = 1, . . . , r. All couplings in the Lagrangian (2.3) 
depend on the scalar fields and are determined in terms of the prepotential F(4>). The 
kinetic term of the scalars is given by a non-linear sigma model with metric 

d 2 F 

G AB (<f), <f>) = Im [tab {4>)} tab = d A d B F = q^aq^b ■ ( 2 - 4 ) 

A metric determined in such a way in terms of a holomorphic prepotential is called (rigid) 
special Kahler [37]. Its Kahler potential is of the special form 

K = Im [d A F4> A ] . (2.5) 

The metric Gab also appears as a scalar dependent gauge coupling in front of the 
kinetic term of the U(l) gaugefields in (2.3). Furthermore there is the term involving the 
dual fields *F A , note that in our notation 

* F A»v = ^ upa F A . (2.6) 
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This term couples to the scalars through 



®ab = Re [tab] ■ (2.7) 

The theory thus also contains a matrix of scalar dependent ^-angles &ab- In summary, 
the matrix tab = ©ab + iGab plays the role of the complexified gauge couplings. 

2.3 Duality structure 

The electromagnetic duality of Maxwell's equations generalizes to theories of the form 
(2.3), where it is part of a larger duality group [38]. To make this duality manifest it is 
convenient to formally extend the U(l) field strengths F A to a doublet 

F a = , Fa = Re [t ab {F b + i * F B )] . (2.8) 

Here the index a = 1, . . . , 2r runs over both upper and lower indices A = 1, . . . , r. To create 
the doublet we introduced a new component Fa, that is defined in such a way that the 
Bianchi identity together with the equations of motion following from (2.3) take a simple 
and elegant form in doublet notation: 

dF a = . (2.9) 

This form of the equations is manifestly invariant under GL(2r, M.) transformations. One 
should however keep in mind that F a can't take values in all of the abstract 2r dimensional 
vector space, since there is a relation between the top and bottom components. We can 
express this relation (2.8) in a covariant way as 

*F«=T\F» with X\=\ GAB r t[f^A Q]AB ) ■ (2-10) 




B 



One can check that I is a complex structure, I = —1, so the above constraint is a 
generalized imaginary self-duality condition. Invariance under duality thus also imposes 
that this constraint remains invariant, i.e. *F' a = I' a t,F' b . Some algebra reveals that this 
is only the case for transformations in the subgroup Sp(2r, M)cGL(2r, M) [38]. 

Note that although the field strengths F a transform in the fundamental representation, 
the couplings tab transform under the natural action of Sp(2r, K) on the Siegel upper half- 
space. Split the symplectic matrix M G Sp(2r, M) into doublet notation as 



M = „ with A T C - C T A = , B T D - D T B = and A T D - C T B = 1 . 




Then the couplings transform as follows under application of M: 

t' ab =[(At + B)(Ct + D)- 1 ] ab . (2.11) 
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These couplings are however not independent fields but determined in terms of the scalar 
fields as the second derivative of the prepotential: tab = OaObJ 7 . To make their transfor- 
mation consistent with this definition one introduces a dual scalar 

<\>A = dAT- (2.12) 

One can now check that the transformation (2.11) follows if we assume that the scalars <j) A 
transform together with the <j>A in a fundamental doublet of the Sp(2r, M) duality group: 




(2.13) 



Because duality mixes the scalar fields with duals defined in terms of the prepotential T it 
follows that this potential itself must transform non-trivially as well. The transformation 
can be obtained by integration from (2.13), the result is 

F' = F+ 1 -4>b [C T A] BC <p c + <p B [C t B] b c4>° + \<f> B [D T B] BC <f ■ (2.14) 
This can be written in the slightly more manageable form as 

r-\M A = T- \<Pa<P A . (2.15) 

As these formulae make clear, and as is important to note, the action (2.3) is in general 
not invariant under duality transformations. It are only the equations of motion following 
from the action that are invariant. As we will discuss in section 4 there exists another 
variation principle leading to the same equations of motion that is more manifestly duality 
covariant. 



3. BPS solutions 



In this section we derive the BPS equations from the Hamiltonian in the standard way and 
show how they are solved in terms of 2r real harmonic functions. The connection between 
these harmonic functions and the r complex scalars is made through a Legendre transform 
on the imaginary part of the prepotential, called the Hesse potential. 

3.1 Hamiltonian formalism 



As usual in a Hamiltonian treatment we need to split the coordinates in time and space, 
i.e. x M = (t, x). Und 
and magnetic compi 
from (2.3) are then 



i.e. x^ = (t, x). Under this split we can then decompose the field strength F A V into electric 
and magnetic components E A and B A . The momenta conjugate to <f> A and A A computed 



7TA = = 2~ Im ^ AB ^ 9t< ^ B ' ^ 3 ' 1 ^ 

Sa = T^lTl = Im Wb{E b + %B B )\ . (3.2) 
5d f A A L J 
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We use the notation Ba for the momentum conjugate to A A , as one can check that (3.2) 
is equal to the magnetic component of Fa, as defined in (2.8). 

Using these ingredients one computes the Hamiltonian corresponding to (2.3): 

H = J d 3 x hm [Gab] (e a ■ E B + B A ■ B B + d t <j) A d t 4> B + V<f> A ■ V0 B ) . (3.3) 

Remember that, as in any U(l) gauge theory, the equations of motion following from this 
Hamiltonian need to be supplemented by the Bianchi identity and the Gauss constraint: 

V-B A = 0, V • 5 A = . (3.4) 

Remark that for the moment we ignore possible sources in these equations as we will focus 
on the solutions in vacuum regions of space-time. We comment shortly on the introduction 
of sources in section 4.3. 

3.2 BPS equations 

To find the BPS equations we rewrite the Hamiltonian (3.3) as a sum of squares and a 
total derivative: 

H = J d 3 x^\E + iB- e - ia V(f)\ 2 + -\d t (f)\ 2 + V-Re e~ ia {(t> A B A - 4> A B A ) ■ (3.5) 

Here we used the norm \X\ 2 = GabX a ■ X B , and introduced an arbitrary constant phase 
a. The total derivative term gives us the central charge 

A j5 j. 3A\ ±A„ j.oo A 



Z = / dn- ^ A B A - <PaB a ) = (p^q A ~ <t>A°p A ■ (3.6) 

In the last step we assumed the scalar fields to asymptote to constant values (j)^ at large 
radius and introduced 'electric' and magnetic charges 



QA 



dn-B A , P A = dn-B A . (3.7) 



Due to the positivity of all terms in (3.5) and the freedom to choose a, the BPS bound 
follows directly: 

H>\Z\. (3.8) 

Solutions that preserve half of the supersymmetry saturate the bound (3.8) and are 
known as BPS solutions. Instead of looking at the Killing spinor equations we can also 
directly check what the conditions are that the Hamiltonian equals the norm of the central 
charge. From (3.5) it follows that this is the case if and only if the following BPS equations 
are satisfied: 

d t cp A = , (3.9) 
arg Z = a , (3.10) 
E A + iB A - e~ ia V(p A = . (3.11) 

The first equation tells us that all BPS solutions should be stationary and the second 
equation fixes the constant a that appears in the last. This last BPS equation (3.11) is 
a set of r first order complex equations. It can be combined with the Gauss and Bianchi 
equations into 2r real Laplace equations, as we will now show. 
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3.3 Harmonic conditions 

To start we split (3.11) in its real and imaginary parts 

Vim [e- ia (j) A ] = B A , (3.12) 
VRe[e" iQ ^] = E A . (3.13) 

The first equation (3.12) is readily solved, as together with the the Bianchi identity (3.4) 
it implies that the imaginary parts of the scalars are harmonic functions: 

Aim [e- ia 4> A ] = 0. (3.14) 

The second equation (3.13) can't be treated similarly since in general E A doesn't need to 
be diver genceless. The way to proceed is to multiply the BPS equation (3.11) by tab'- 

t ab {E b + iB B - e- lQ V^) = t ab (E b + iB B ) - vV^a) = . (3.15) 

This equation now has as imaginary part 

Vim [e- ia (j)A] = B A . (3.16) 

Since Ba is indeed divergenceless by the Gauss constraint (3.4) we find that also the 
imaginary parts of the dual scalars are harmonic: 

Aim [e~ ia (t)A] =0. (3.17) 

So the BPS equations imply 2r real Laplace equations (3.14, 3.17) on R 3 . Given a set of 
appropriate 2 sources these can be solved through standard techniques. The electromagnetic 
fields are then simply the gradients of these harmonic functions. What still has to be 
clarified however is how the r complex scalars are related to the 2r real solutions to the 
Laplace equation. In the next subsection we show that this relation is naturally interpreted 
as a Legendre transform. 

3.4 Legendre transform 

The motivation for the discussion in this subsection comes form the analysis of the BPS 
equations made above. It is however important to point out that all formulae in this 
subsection are valid off-shell. In section 4 we will use this to reformulate the theory directly 
on the level of the action. 

The question we want to address is if and how we can express the r complex scalars 
(j) A as a function of the 2r real scalars 3 Im [e~ ta (/>^] and Im [e -10 ^] . To lighten notation 
and clarify the discussion we will explicitly split up the complex scalars into their real and 
imaginary parts: 

e~ ia 4> A = X A + iY A , (3.18) 
e- ia <t> A = X A + iY A . (3.19) 

2 As we discuss in section 4.3, the sources can interact with one another, which will put constraints on 
what sources are allowed for a BPS solution. 

3 In principle we can take a to be an arbitrary phase. If we want to relate this formalism to the BPS 
solutions we will of course make the identification (3.10) 
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From this split it is clear that the problem to find <f) A (Y A , Ya) reduces to finding X A (Y A , Ya)- 
To start remember that 4>A is defined in terms of cj) A through the prepotential i.e. 
(j)A = Furthermore the prepotential is holomorphic, which implies it satisfies the 

Cauchy-Riemann conditions 

dT_ _ ia dT_ _ ia dT_ 

d<j> A ~ e dx A ~ ie dY A • [6 - 2U) 

Combining these two facts one finds the relations between the different real scalars: 

dim \e~ 2ia T\ 
aim \e- 2ia T] 

The form of these relations shows that we can think of (Y A ,X A ) and {X A ,Y A ) as pairs 
of conjugate variables and we can trade one component of a pair for the other one by 
performing a Legendre transform. By our motivation from the BPS equations we will 
prefer a formulation where the Y's are the 2r real variables, so we define the Legendre 
transform 

S(Y A , Y B ) = Y C X C - Im [e- 2ia F{X A , Y A )] . (3.23) 

In this transform the scalars X A are now functions of the 2r real scalars (Y a ,Ya), found 
by inverting (3.22). If we consider S as a given, then we can find those functions simply 
by taking a derivative: 

*' 4 = aSi £ M < 3 - 24 > 

This formula provides, at least formally, the solution to our problem, we can now write 

4>\Y^ Y C ) = e <° [ 8S{Y d yJ C) + iY A ^J . (3.25) 

This expression (3.25) is exactly what we need to complete our solution of the BPS 
equations, since these simply imply that the Y's are all harmonic (3.14, 3.17) . So given the 
appropriate harmonic functions one can plug them into (3.25) to find the corresponding 
complex scalars. 

Finally let us point out that one can also directly compute the dual scalars in terms 
of the Y variables: 

dS(Y B ,Y c ) ( dS(Y A ,Y A ) \ ( . 

X A = 4>A = e ^ — A + iY A ) . (3.26) 

The minus sign that appears here is important and is directly related to the symplectic 
structure on moduli space, as we will discuss in the next section. 
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4. The Hesse potential: duality and special geometry 

In our study of the BPS equations in the previous section we were naturally led to introduce 
a function S as the Legendre transform of the imaginary part of the prepotential J-, see 
(3.23). In this section we will show that the Hesse potential S is invariant under duality 
and that the special geometry of the theory can be elegantly expressed in terms of it. Most 
of this will be review of the discussions in [30, 31, 34]. In the second part of this section 
we show how the theory can be expressed in a duality covariant description that seems 
most naturally connected to the BPS equations and solutions and we comment on the 
constraints on sources. 



4.1 The Hesse potential 

As we will show in this subsection, S plays a role in special geometry that is very analogous 
to that of the prepotential J 7 , the difference being that J- makes the complex structure of 
the manifold manifest while a description in terms of S makes the symplectic and duality 
structure manifest. This will present itself in two ways, first of all we will show that S is 
duality invariant and secondly that in Darboux coordinates the metric is simply given by 
the Hessian of S. Due to this last property S has been called the Hesse potential [27, 29]. 

To check the duality invariance of S only takes a few lines of algebra. It follows from 
combining the definition (3.23) in terms of the prepotential J 7 , with the transformation 
properties (2.13, 2.14): 

S' = Y' A X' A - Im [e~ 2ia T'] (4.1) 
= S + \ (X A Y A - X A Y A ) - \ (X' A Y ,A - X' A Y' A ) (4.2) 
= S. (4.3) 

The last step follows because X A Y A — X A Y A = X' A Y' A — X' A Y' A is an invariant under the 
symplectic duality group Sp(2r, R). This is most clear if we introduce doublet notation 
and a symplectic form: 



* a = Z A A > . n= " ■ ( 4 - 4 ) 






In this notation X A Y A — X A Y A = £l a bY a X b which is manifestly invariant under symplectic 
transformations. Note that this invariant is nothing but the Kahler potential (2.5), i.e. 

K = n ab x a Y b . 

The invariance of S under duality rotations suggests that it is closely related to the 
symplectic structure of the theory. Indeed, as we will now describe it plays a role similar 
to the prepotential, in that its second derivatives give the metric on the scalar manifold. 
However, since S is a function of the 2r real coordinates Y a instead of the r complex cj) A , 
the matrix of second derivatives is nothing but the Hessian matrix of S. 
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To make this explicit one can start from the scalar metric in complex coordinates (2.4) 
and rewrite it using the relations (3.24, 3.26): 



G A Bd4> A d(t) B = Im [d<f> A d<j> A ] 

= dY A dX A - dY A dX A (4.5) 

-dY a dY b . 



8 2 S 



Qyagyb 

In summary the scalar metric in the 2r real coordinates Y a is indeed the Hessian matrix 
of S, we will use the following notation: 



d 2 s 



Jab 



8Y a dY b ' 



(4.6) 



What is interesting is that the metric Q ao is actually a rather familiar object. Remember 
that we have both a symplectic form fi a & and a natural complex structure I a b defined in 
(2.10), that are furthermore compatible, i.e. Q a b = £l c } } T c a Z d i ) . This implies that together 
they define a metric: 

(q—i —G~ 1 \ 

-GG^ G + BG-1&) ■ (47) 

As we already implied in our notation, this metric is nothing but the Hessian metric (4.6). 
This was shown in e.g. [32, 39, 34], but for completeness and since we will use this result 
later on in the paper, we have added a derivation of this fact using our notation in appendix 
A. 

4.2 Duality covariant formalism 

In the discussion of the BPS equations and their solutions in section 3 we found that 
those are most naturally expressed in terms of the 2r real scalars Y a and the magnetic 
fields B a . To realize this we however had to first perform some algebra and derive the 
function S from the prepotential T . Using the results of the previous subsection we will 
now present an alternative formulation of the theory where these variables appear as the 
fundamental degrees of freedom in the action and the simplest form of the BPS equations 
follows immediately. By construction this new action will be manifestly covariant under 
duality. In appendix B we show in detail that the new formulation leads to the same 
equations of motion 4 as the more standard description in terms of the action (2.3). 

In the original description (2.3) there are r complex scalars <p A and r 4d gaugefields 
A A as fundamental variables. The gauge potentials only enter the action through the field 
strengths F A = dA A , which is thus manifestly invariant under U(l) r gauge transformations 
A A — > A A + d^A. In this formulation duality transformations act in a nontrivial non- 
local way. For such Abelian gauge theories there exists a nice alternative formulation that 



4 Although in this paper we are mainly interested in BPS solutions and equations, the equivalence of the 
theories is checked for the complete set of second order equations of motion without any further assumptions 
on BPS or other conditions. 
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makes duality manifestly covariant and furthermore implements it as a local transformation 
[40, 41]. Combined with the field redefinition (4.5) in the scalar sector this formalism is 
naturally extended to any N = 2 low energy gauge theory. 

In this alternative description we will take as fundamental variables 2r real scalars Y a , 
2r 3d vectors A a and 2r functions Aq. Note the naive 'doubling of degrees of freedom' in 
the electromagnetic sector with respect to the usual formulation (2.3). Using the vector 
potentials we define the following magnetic and electric fields: 

B a = V x A a , E a = d t A a - VA a . (4.8) 

By construction these satisfy the following 'Bianchi-identities': 

V • B a = , d t B a = V x E a . (4.9) 

Given the symplectic form £l ab = —e ab an d a scalar dependent metric Q a b(Y), we write 
down the following Lagrangian: 

l = - x - (n ab E a ■ B b + g ab B a ■ B b + g ab vY a ■ vY b - g ab d t Y a d t Y b ) . (4. 10) 

First of all it is interesting to note that this Lagrangian has the following two gauge 
invariances 

A a ^A a + A a E a ^E a -VA a , L^L + ^V-(U ab A^B b ), (4.11) 

A a ^A a + VA a => E a ^E a + d t VA a , L —?■ L — • (U ab d t A a B b ) . (4.12) 

Remark that the electric fields in this formulation are not gauge invariant, contrary to 
what we are used to in Maxwell's formulation of electromagnetism. Furthermore the first 
line shows us that the fields Aq are pure gauge and can be put to an arbitrary value. This 
is consistent with the fact that they appear as a total derivative in the Lagrangian (4.10), 
hence varying with respect to them will not lead to any non-trivial field equation. Again 
this is different from the standard Maxwell action. 

Furthermore under the Sp(2r, M) duality group all fields now simply transform in the 
vector representation, i.e. for M a b £ Sp(2r, M) we have f a — > Ad a b f b . Note that this is 
also true for the vector potentials A a , so these now transform in a manifestly local way. 
The Lagrangian (4.10) is in general not invariant under these duality transformations, only 
for special choices of Q ab this will be the case. However, it immediately follows from the 
covariance of all fields in the action that the set of equations of motion is invariant under 
duality transformations. 

We deliberately used the same notation for the fields appearing in the Lagrangian 
(4.10) as for those we introduced earlier in the text. As is shown in appendix B, under 
these identifications the equations of motion following from (4.10) are equivalent to those 
of the action (2.3), and so both Lagrangians are nothing but different descriptions of the 
same physical theory. 
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One can derive the BPS equations directly in this new description. To do so let us 
again go to the Hamiltonian formulation, the momenta conjugate to the new variables are 

n a = -^r- = -ln ab B b , n a = g ab d t Y b . (4. 13) 

5d t A a 2 

The Hamiltonian density is then computed to be 

U = X -Q ah (B a ■ B b + VY a ■ VY b + d t Y a d t Y b ^j . (4.14) 

Note that contrary to the Lagrangians the two Hamiltonians (3.3, 4.14) are actually equal 
on-shell, as expected from their physical interpretation as the energy density of the solu- 
tion. We can also rewrite the Hamiltonian density (4.14) as a sum of squares plus a total 
derivative 

ft = \Sab(B a - VY a )(B b - VY b ) + \p ah dtY a d t Y b + V • (d a S B a ) . (4.15) 
This expression immediately implies the following bound: 

H> [ d a SB a -dn. (4.16) 

In the next section we will discuss how this bound is nothing but the standard BPS bound 
in our new notation. Demanding the bound to be saturated gives the BPS equations, that 
now take a very simple form: 

d t Y a = 0, B a = VY a , AY a = . (4.17) 

Here the third equation follows from the second by the Bianchi identities (4.9). Of course 
the equations (4.17) are equivalent to the BPS equations derived in section 3. 

4.3 A comment on sources 

Since the BPS solutions are simply determined in terms of harmonic functions one might 
naively assume that one can take linear superpositions of any number of solutions. This 
is however not the case when these carry mutually non-local charges. The subtlety comes 
from the precise treatment of the sources to these harmonic functions. We leave a general 
and detailed discussion to future work, but illustrate the issue in the case of static pointlike 
sources. In this case we need to modify the equations (4.17) to AY a = ^2i^i^ 3 (x — Xi), 
which gives 

Ya = EA\+ Y ~- ( 4 - 18 ) 

i 

To fully take this effect into account we also need to add the energetic contribution of the 
sources to the Hamiltonian (4.14), i.e. 

H tot = H + Y,\ z i\ 63 (x-^> Z i = n ab a (x)T b . (4.19) 

i 
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When we rewrite T~L as a total square as in (4.15), the use of the Bianchi identities now give 
rises to an additional term — ^ Re [e ia Zi\ 5 3 (x — x~i) . To have a solution that saturates 
the BPS bound this term needs to cancel against the source contribution in (4.19), which 
leads to constraints on the allowed relative positions of the sources: 

E^=W!. (4.20) 

3 J 

These are the well known stability conditions [18] that lead to the non-trivial bound state 
structure of generic J\f = 2 BPS solutions. Note that apart from this algebraic constraint on 
possible source locations our analysis before was completely general and is valid everywhere 
in the vacuum away from the sources. In the rest of the paper we focus again on this vacuum 
behavior of the solutions. 

5. Generalized attractor equations and Hesse flow 

BPS solutions to low energy N = 2 supergravity and gauge theories have been especially 
well studied in the spherically symmetric case. In this case there is only dependence on 
a single radial coordinate, and it turns out that the BPS equations imply a simple radial 
flow for a local version of the central charge. This phenomenon is probably best known in 
the supergravity context, where it goes under the name of the attractor mechanism. This 
result gives a nice physical interpretation to the BPS solution. We start out with boundary 
conditions that set the central charge at infinity, and as we then move radially inward this 
central charge becomes smaller and smaller, reaching its minimum at a certain radius r*. 

In this subsection we want to show how this formalism generalizes to the cases without 
spherical symmetry, where in general there is dependence on all three spatial coordinates. 
As we will argue, the standard notion of local central charge will need to be generalized to 
a central charge vector field. We we will then show that the BPS equations imply that the 
real part of this central charge vector field becomes the gradient of a potential, and that 
this potential is none other than the Hesse potential S evaluated on the solution. Moreover, 
we will see that this equation is only one in a set of four, that take the form of Maxwell's 
equations. The role of electric charge will be played by the energy density in the solution, 
while there is also a magnetic component to the central charge vector field that couples to 
the momentum current in the solution. 

5.1 The central charge vector field 

So far we have defined the central charge of a solution (3.6) as the complex number 

Z = P A <f>A-QA<l>i. (5.1) 

As is well known [42], this is the quantity showing up in the superalgebra of the theory. 
More generally we can define a bilinear central charge function that maps any duality 
covariant real vector V a to a complex number for each value of the scalar moduli cf> A , i.e. 

z : R 2r x M -»• C : (V,<f>) z(V, (j)) = tt ab (j) a V b , (5.2) 
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where we again used the notation <p a = (<fi ,8aJ~) as in (2.13). So in our terminology the 
central charge of a solution is nothing but the image under the central charge function, of 
the total charge T a = (qA,P A ) and the asymptotic values of the scalars in that solution, i.e. 
Z = z(T, cpoo). In previous studies of BPS solutions it became clear that a closely related 
notion, that of a local central charge, is useful in understanding the physics of the solution. 
Given a scalar field configuration, i.e. (p A {x) and a total charge r a , one defines 

Z(x) = z(T, 4>{x)) = n ab (P a (x)T b . (5.3) 

By definition the local central charge evaluated at infinity equals the global central charge, 
i.e. limA-5.00 Z(Xx) = Z. In the spherically symmetric case the local central charge only 
depends on the radial coordinate, i.e. Z = Z(r). It is then a well known result that 
the BPS equations have a simple and very physical interpretation in terms of this central 
charge. In the gauge theory setting the BPS equations imply [21, 18]: 

Im [e ia Z] = and d r \Z\ = ^G AB d A Zd B Z > . (5.4) 

Remember that a = arg Z, so the first equation implies the local central charge has a 
position independent constant phase. The second equation then tells us that the norm of 
the central charge decreases with respect to the radial coordinate. A more precise analysis 
reveals that the minimum is reached at finite radius. What is important is that these 
equations only hold in the spherically symmetric case. As far as we are aware there exist 
no simple to interpret equations for the local central charge when there is no spherical 
symmetry and non-mutually local charges are present. 

We will now propose a more general notion of 'local central charge', for which a nice 
physical interpretation of the BPS equations will appear for all BPS solutions and that in 
the spherically symmetric case reduces to (5.4). 

When looking to generalize Z = z(T,cj)(x)) the first thing that comes to mind is 
to also replace the first argument by a space-time dependent quantity. A simple guess 
would be to replace r a — > Y a {x), this is interesting as this new 'local central charge', 
Z(x) = z(Y(x), 4>(x)) = Q. a b<p a {x)Y b {x) has a fixed phase by construction: 

= Q ab Y a Y b = . (5.5) 

So Z seems to at least generalize the first equation of (5.4) to generic solutions. However 
there seems to be no nice equation describing the spatial dependence of the norm of Z, 
which would generalize the second equation in (5.4). Furthermore there is the additional 
problem that in the special case with spherically symmetry it becomes clear that Z is a 
physically different object than Z. Spherical symmetry implies that Y a = + and 
so Z = — + z(Y OQ , 4>(x)) in that case. The second term is non-trivial and so Z satisfies 
different equations than Z, even in the simple spherically symmetric case. Hence this is 
not the generalization we are looking for. 

A more fruitful proposal is to replace Z by a 3- vector field Z that we define as follows: 



Im 



Z(x) = z{B(x),(P{x)) = Q ab (P a (x)B b (x) 



(5.6) 



Note that here B a is the symplectic vector containing both the electric and magnetic fields 
(2.8). A first good property of this central charge vector field? Z, that justifies to call it 
a generalization of Z, is that in the spherically symmetric case both are essentially the 
same. Assuming spherically symmetry the electromagnetic fields have the form B a = 
and thus Z = -gr. So in this special case Z shares the same constant phase with Z and 
the norms of the two quantities are proportional by a monotonous radial factor and thus 
exhibit the same attractor behavior. 

In the more general non-spherically symmetric case, especially when multiple non- 
mutually local charges are present, Z is no longer simply related to Z. However, as we 
will show explicitly in the next subsection, it is Z that satisfies simple first order equations 
in any context, with no additional symmetry assumptions. This suggests Z is really the 
natural local generalization of the global central charge Z defined at infinity. 

5.2 Hesse flow and central magnetism 

Let us now show how one can generalize the equations (5.4) to cases without spherical 
symmetry and non- mutually- local charges, making use of the central charge vector field Z 
defined in (5.6). Even in this more general case, if we assume all charges to be contained in 
some finite region, then near infinity the field strengths will asymptotically approach those 
of a spherically symmetric distribution. So it follows that the phase of the central charge 
vector field at infinity is equal to that of the global central charge (3.6). In formulas 



lim Im 

A— >oo 



e - ia Z(\x) =0 with a = arg(Z). (5.7) 



But once we move away from the boundary at infinity 'the' phase of Z will become position 
dependent. To be more precise, each of the three components actually develops an inde- 
pendent position dependent phase. So the description in terms of the norm and the phase, 
that we used to formulate (5.4), seems not so useful for generic BPS solutions. Rather it 
will prove natural to split up the complex central charge vector Z in its real and imaginary 

components, Re[e -JQ iT] and Im[e~ lOL Z\. Note that this split becomes related to the norm 

~* \z\ 

and phase when the solution has spherical symmetry, then one finds Re[e~* a i?] = and 
lm[e- ia Z] = 0. 

Let us start our analysis with the real part. Simply writing out its definition from (5.6) 
in terms of the symplectic coordinates Y a and their Hesse potential S that we introduced 
in the previous sections, one finds Re[e~* a 2*] = B a d a S. In this form it becomes a simple 
observation that the BPS equations B a = VY a imply a gradient flow equation: 



Rc 



VS. (5i 



Note that here the Hesse potential S is to be considered as a function on R 3 , evaluated as 
S(Y(x)). Equation (5.8) shows that the Hesse potential, apart from being an ingredient in 



Note that this vectorfield is analogous to the graviphoton field strength in supergravity, as it is essentially 
the unique duality invariant contraction of the scalars and electro-magnetic fields. We thank T. Mohaupt 
for pointing this out to us. 



- 16 - 



the special geometry of the theory, also gets a physical interpretation on BPS solutions, as a 
potential for the real part of the central charge vector field. This gradient flow in 5, which 
we will refer to as Hesse flow, becomes the well know attractor flow in the spherically 
symmetric case, but is more general and as we just showed holds for all BPS solutions 
irrespective of symmetries or charges. 

The Hesse flow equation (5.8) is however only the first manifestation of a bigger un- 
derlying structure characterizing BPS solutions. A rather trivial consequence is that for 
BPS solutions the real part of the central charge vector field is irrotational, i.e. 



V x Re 



0. 



(5.9) 



More interesting is to compute its divergence. Evaluating it for vacuum BPS solutions one 
finds 



V • Re 



n. 



(5.10) 



Here T~L is the energy density (4.14) of the BPS solution. Note that this equation is not at 
all unexpected, it is simply the local version of the well known BPS condition H = \Z\, 
and it is an immediate consequence of the expression (4.15). Combining it with the Hesse 
flow equation (5.8) reveals an intriguing relation however: 



As = n. 



(5.11) 



Here we see that the Hesse potential evaluated on a BPS solution satisfies Poisson's equation 
sourced by the energy density. In physical terms, for a BPS solution the Hesse potential is 
non other than the gravitational potential! Furthermore (5.8) then gives the real part of 
the central charge vector field the interpretation of the gravitational field. 

After uncovering these rather elegant equations for the real part of the central charge 
vector field, let us analyse the imaginary part. Although this is trivially zero in spherically 
symmetric solutions, that is not the case in general. Again one starts from the definition 
(5.6) and rewrites things in terms of the real scalars Y a : Im[e~' ia 2] = Q a bY a B b . Since the 
BPS equations read B a = VY a and the Y a are harmonic it follows immediately that 



V • Im 



0. 



(5.12) 



Furthermore the BPS equations imply an elegant expression for the curl as well. One 
computes 



V x Im 



2V. 



(5.13) 



Here V = ^ft a bB a x B a is nothing but the Poynting vector field of the solution. 

A very nice structure unfolds when we put all these expressions together. The calcu- 
lations above show that the BPS equations imply that the central charge vector field (5.6) 
satisfies the following equations: 



V-Re 



-Z 



He 



V x Re 



V • Im 

V x Im 



(5.14) 
(5.15) 
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These are none other than Maxwell's equations, with (half of) the real part of the central 
charge vector field playing the role of 'electric field', while the imaginary part behaves 
magnetically. Note that the role of 'electric charge density' is played by the energy density 
of the electromagnetic fields in the theory, Ji e . m . = for BPS solutions. The current 
V that sources the 'magnetic' imaginary part is none other than the momentum current. 
Indeed these equations are Lorentz invariant as H e .m. an d V nicely combine in the energy- 
momentum 4-vector. 

It is straightforward to check that in the spherically symmetric case the magnetic 
equation becomes trivial, while the electric flow equation becomes the attractor equation 
(5.4) (up to a trivial radial rescaling). 
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A. Duality covariant metric and the Hesse potential 

In this appendix we show that the duality covariant metric introduced in (4.7) is the Hessian 
of the potential S defined in (3.23), i.e. Q a b = Qylgyb ? with S = YaX a — Im \e~ 2ia T\ and 

e=( °~\ - G "' e , V (A.i) 



Before we start the derivation of this fact, let us for convenience first summarize some 
definitions and identities from the main body of the paper. Remember that X A = -^x, 
Xa = ~~§^ an d we use the notation tab = qJa^b — @AB + iGab- Furthermore a useful 
form of the Cauchy-Riemann equations for any holomorphic function H((p), (J) = e ia (X+iY) 
is 

dH_ _ ia dH_ _ . jg 9H 

dcjA ~ 6 dx A ~ ie 8Y A ' [ > 

Combining the different identities above one finds that 

9Ya 9Xa A (AV 
W^ = dx^ = ° AB ' (A - 3) 

QYa 9Xa r ( a a\ 

dxB = W^ = GAB - (A - 4) 

Let us furthermore consider the coordinate transformation (X A ,Y A ) — > (Ya,Y a ), the 



Jacobian of this transformation and it's inverse are given by: 

/ 8Y A dY A \ ( dx A dx A \ 

t_ dX B 8Y B 7-1 _ dY B 8Y B CA K\ 

- [ S A ) ' - ^ 5 A ) ■ 
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Because of the fact that J J =1 the following relations follow: 



dY A dX B c 

dx*Wc =SA > (A - 6) 

8Y A dX B dY A 



dX B dY c " dY c 



(A.7) 



We can now use the results above to compute the Hessian of S in terms of the original 
couplings r. It is a matter of simple algebra to check that 

Computing the remaining components of the Hessian is a little more subtle, one needs to 
take into account that 

d 2 s dx A ax A ax A dY c 



dY A dY B dY B dY B dY c dY B ' 
We can now again apply the various identities above to find: 



f)2c 

WAW s=G AB +[QG~'e] AB . (A.10) 
Summarized in matrix notation the results look as follows 

Hpqs rc\_ dY A dY B dY A 8YB _ 1* -tr V C B ail) 

{ ) " l^fe wMy* J " \-^cG CB g ab + e^^e M J • (A - n) 

By comparing with (A.l) we see that indeed, as we set out to show, Hess(5) = Q. Using 
the index notation Y a = (Y A ,Y A ) this can be written as Q a f, = gylgyb ■ 



B. Equivalence of duality covariant description 

In the main text we introduced two Lagrangians, (2.3) and (4.10), (and their correspond- 
ing Hamiltonians, (3.3) and (4.14)) to describe a single low energy J\f = 2 gauge theory. 
The Lagrangians are functions of different variables and are extremized with respect to 
two different variational principles. But as we will explicitly show in this appendix, they 
lead to equivalent equations of motion if one identifies their fundamental fields appropri- 
ately. For the electromagnetic sector of our theories both the formulation of the duality 
covariant Lagrangian and the check on its equations of motion are a rather straightforward 
generalization of the formalism of [40, 41]. In our particular case there is however also 
a scalar sector that couples non-trivially to the electromagnetic fields and we show that 
even including those, the complete sets of equations of motion of the two Lagrangians are 
equivalent. 
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B.l The two descriptions 

Let us for clearness start by comparing the two Lagrangians and their fundamental vari- 
ables, at the same time splitting them up into pieces that play different roles. 

The first, most standard and manifestly Lorentz invariant, Lagrangian (2.3) is 

Lx{^Ai) = Lr^.^ + if^^), (B.l) 



= -±G AB F A F B ^ - \&abF* *F b »» , (B.2) 

Lf kin - = -\G AB d^ A d^ B . (B.3) 

Its fundamental variables are r complex scalars <fi A and r 4d gaugefields A A . We will call 
this theory 'description 1'. 

The second, manifestly duality covariant, Lagrangian (4.10) is 

L 2 (Y a ,A a ,A%) = L c 2 m iA a ,A%) + L™- +sc -{A a ,Y a ) + Lf- kin -(Y a ), (B.4) 
L c 2 m - = -hl ah E a -B\ (B.5) 
L m.+sc. = _^ gab B«.B\ (B.6) 

^ S 2 C ' kin - = \Qabd t Y a d t Y b - l -G ab VY a ■ VY b . (B.7) 

Here the fundamental variables are 2r real scalars Y a , 2r 3d vector potentials A a and 2r 
potentials Aq. Note that these last are pure gauge and their variation doesn't lead to a 
non-trivial field equation, see section 4 for some details. This second theory we will refer 
to as 'description 2'. 

We will now show that the field equations of these two theories are exactly the same 
and the two are different descriptions of the same physics, provided we make the following 
identifications: 

Y a = 



B a = \ %Ca I with Ba = G AB E B + Q AB B B , (B. 



lab 




-[G-i@) A B \ 

u G AB +[e G -i@] AB ) ■ 

It is important to note that for each field and coupling in description 1, we have now 
identified a corresponding physical field or coupling in description 2. There remains however 
a field in description 2, i.e. E a , that so far (i.e. off-shell) has no counterpart in description 
1. We cannot link it to a physical field in description 1 without introducing a constraint, 
since above we have already associated a field of description 2 to each field of description 
1. As we will discuss, the equations of motion following from the respective Lagrangians 
provide exactly this constraint and so on-shell the two theories are equivalent. Furthermore, 
as we will stress and explain below, it is important that the identifications (B.8) are made 
only after the respective Lagrangians have been varied. 
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B.2 Electromagnetic field equations 

Let us start by showing that the field equations for the electromagnetic fields are equivalent. 

In description 1 the independent fundamental electromagnetic fields are E A and B A . 
As we described in the main text in section 2.3, the field equations of description 1 take 
the following form: 

• Bianchi identities 

dF A = <^ V • B A = , d t B A = V x E A . (B.9) 

• Definitions 

B = 1E B A = G AB E B + e AB B B , E A = -G AB B B + @ A bE b . (B.10) 

• Equations of motion 

dF A = & V-B A = 0, d t B A = VxE A . (B.ll) 

In description 2 there are 2r independent magnetic fields B a , furthermore there are as 
many electric fields E a . Note that these last are however not gauge invariant and so not 
all of their degrees of freedom are physical. The field equations for these electromagnetic 
fields are the Bianchi identities (4.9) together with the equations of motion following from 
varying the Lagrangian (B.4) with respect to A a : 

• Bianchi identities 

V-B a = 0, d t B a = VxE a . (B.12) 

• Equations of motion 

V x (n ab E b + g ab B b ^j = . (B.13) 

Let us now show that the equations (B.12, B.13) are equivalent to (B.9, B.10, B.ll) 
assuming the identifications (B.8). Using the relation between the symplectic form f2, the 
metric Q and the complex structure X given in (4.7) we can solve equation (B.13) as 

B a = l a b E b + W 6 . (B.14) 

Due to the gauge freedom (4.11) we can however always put "$> b = 0. Using the identification 
(B.8) that B a = (Ba, B a ), equation (B.13) then implies together with the definition (B.10) 
that E a = (Ea,E a ). It then immediately follows that equations (B.12) are equivalent to 
equations (B.9, B.ll). 
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B.3 Scalar field equations 

To compare the scalar field equations of the two theories (B.l, B.4) it is best to split the 
problem into two parts. Under the identification (B.8) we can think of the complex scalars 
as a function of the real ones: <ft A = cp A (Y a ). As we explicitly checked in section 4 we have 
the following equality: 

Lf- kin -(0 A (Y a )) = Lf- kin -(Y a ) . (B.15) 

As is well known the Euler-Langrange equations remain satisfied under a scalar field redef- 
inition, i.e. 

,5L s 2 c ' kin ' 5Lf- kin - _ f 5Lf- kin - 5Lf- kin -\ d(p A 

9 >*ld^~~8Y^~ {^^d^ ~ ~d^ A ~ J W^ + C - C -- (B ' 16) 

It is important to realize that the same reasoning doesn't work directly for the other 
terms in the Lagrangians. The variation principle for the scalar field equations in descrip- 
tion 1 is to vary the scalars (f> A while holding E A and B A fixed, while in description 2 we 
vary Y a while keeping B a = (B A ,B A ) fixed. Since the relation B A = G AB E B + @ AB B B 
depends explicitly on the scalars we cannot simply relate the parts of the Lagrangians 
involving electromagnetic fields by a field redefinition while preserving the variation prin- 
ciple. However, under such a naive field redefinition using the identifications (B.8) the 
two Lagrangians (B.l, B.4) are actually not equal. We will show below that the difference 
exactly compensates for the changed variation principle and that the two effects nicely 
cancel out, so that 



^m.+sc. / 5L c.m.\ 



5Y a V H A J 9Y a 
This relation then combines with (B.16) to yield 



+ c.c. . (B.17) 



kA 



a &L 2 5L 2 _ SLi SLi\ 

^Sd^Y" 6Y" ~ \ ^5d^ A 5<j> A ) dY- + C ' C - ■ 

Since cf> A {Y a ) is by construction assumed to be invertible we thus find that also the scalar 
equations of motion of the two descriptions are equivalent. 

To show that (B.17) holds, let us start in description 1. Varying L| m ' with respect to 
(j) A while keeping B A and E A fixed one finds 

rre.m. i 

-jjjT = - (d A G BC ) {B B • B c — E B • E c ) - (d A G BC ) E B ■ B c . (B.19) 

Similarly we can compute the variation of L™ +sc ' with respect to Y a while holding B a = 
(B A ,B A ) fixed, the result is simply 

rj-m.+sc. -I 

= ^ (d a G bc ) B b B c . (B.20) 

Now that we have completed the variation principle we are free to use the identifications 
(B.8). After some careful algebra it follows that 

l - (d a g bc ) B b B c = l - (d a G AB ) (b a -b b -e a - e b ) - (d a e AB ) e a -b b . (B.21) 
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To see that (B.17) is indeed correct it is now enough to use the fact that for any real function 
f{</>, <j>) the 'chain rule' reads d a f(cb(Y), 4>{Y)) = ^d a (p A + Mzd a 4> A , and compare (B.20, 
B.21) to (B.19) and it's complex conjugate. 
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